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Abstract 

Continuing the programme of constructing the backreacted solution corresponding to 
smeared anti-D3 branes in the warped deformed conifold, we solve analytically the equations 
governing the space of first-order deformations around this solution. We express the results 
in terms of at most three nested integrals. These are the simplest expressions for the space of 
SU{2) x SU (2) x Z2~invariant deformations, in which the putative solution for smeared anti- 
D3 branes must live. We also explain why one cannot claim to identify this solution without 
fully relating the coefficients of the infrared and ultraviolet expansions of the deformation 
modes. The analytic solution we find is the first step in this direction. 
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1 Introduction and Motivation 



Anti-D3 branes in Klebanov-Strassler (KS) warped deformed conifold throats [1] are a key in- 
gredient in string model building and string cosmology, where they are used both for lifting AdS 
to de Sitter solutions [2], and to construct models of inflation using D3 branes moving in KS-like 
geometries [3]. 

In a previous paper [4], three of the authors solved for the space of SU{2) x SU{2) x Z 2 - 
invariant deformations, in which the would-be solution describing anti-D3 branes smeared at 
the tip of the warped deformed conifold must live. The equations parameterizing this space of 
deformations were found in [5], and by repackaging them and finding their homogeneous solu- 
tions, one can write down the implicit solution to the full space of deformations in terms of 
nested integrals [4]. Performing these nested integrals is no easy task. To analyze the physics of 
the backreacted solution sourced by antibranes and determine whether this solution is dual to a 
metastable vacuum of the dual gauge theory (as conjectured by [6] following a probe analysis), 
the authors of [4] found the explicit infrared and ultraviolet expansions of all the modes 1 . This 
analysis revealed that the candidate anti-D3 brane solution must have a certain infrared singular- 
ity. If this singularity is physical, then the space of solutions implicitly-solved-for in [4] captures 
the first-order backreaction of antibranes in KS. If this singularity is pathological, the analysis 
of [4] implies that the backreaction of anti-D3 branes in KS cannot be treated in perturbation 
theory 2 . 

As discussed in detail in [4] there are many arguments why this singularity may be physical 
or not, but so far there exists no criterion for accepting it that does not force one to also accept 
clearly unphysical singularities like that of the negative-mass Schwarzschild solution 3 . Solving 
the fully-backreacted solution for these anti-D3 branes could settle the issue, but this involves 
solving eight coupled second-order nonlinear differential equations. A more indirect way to 
determine whether this singularity is physical would be to allow it, and then explore whether 
the physics of the resulting solution is consistent with the physics one expects from anti-D3 
branes. To do this, one must set all the other unphysical IR singularities and all the UV non- 
normalizable modes to zero (which was done in [4]), and then relate the various coefficients of 
the infrared and ultraviolet expansions. This will fix all the parameters of the solution in terms 
of one parameter, the number of anti-D3 branes, and yield the relation between the charge 
of the would-be antibranes and the force they exert on a probe D3 brane. Since this force is 
known from KKLMMT [3] and is not screened [11], if the force-to-charge ratio is the correct one 
for anti-D3 branes, this indicates that the singularity is most likely physical, and the solution 
describes indeed the backreaction of anti-D3 branes. If this ratio is not the correct one then the 
solution has nothing to do with antibranes, and is most likely unphysical. Note that obtaining 
this ratio from the solution is a nontrivial process, which involves matching all the infrared and 
ultraviolet integration constants in the expansions of the modes. 

1 The UV expansion of a subset of the invariant modes for the Klebanov-Tseytlin solution has been considered 
in [7] and the IR expansion of a deformation around KS has been looked at in [8] . 

2 A similar result was obtained when investigating the backreaction of putative brane-engineered metastable 
vacua [9], or the backreaction of anti-M2 branes in a warped Stenzel background [10]. 

3 An example of such a criterion is accepting so-called "normalizable," singularities whose energy density 
diverges, but whose energy is finite. 
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One way to obtain this ratio would be to perform numerically all the nested integrals that 
appear in the solution of [4] and relate the corresponding ultraviolet and infrared integration 
constants. However, this gets complicated, in particular because most of the integrals either 
diverge or decay very fast. What simplifies the analysis is the observation that one of the 
modes (£1) is proportional to the warp factor in the Klebanov-Strassler solution. This simplifies 
the expressions of the integrals appearing in the functions £ a , which in the Borokhov-Gubser 
formalism are the "conjugate momenta" to the physical perturbation fields <p a [5]. The second 
key observation is that most of the nested double integrals that appear at various stages of 
solving the Borokhov-Gubser equations can be integrated by parts, and hence can be expressed 
in terms of single integrals. 

In this note we show that by making judicious use of these two observations, one can express all 
eight conjugate momentum functions £ a as linear combinations of simple integrals. Furthermore, 
with the exception of the perturbation to the warp factor, all the other seven perturbation modes 
can be expressed in terms of just two nested integrals. The expression for the perturbed warp 
factor is the most complicated of them all, and can be expressed in terms of three nested integrals. 
This is a huge improvement from the implicit solution of [4], where the most complicated of the 
£ a was expressed in terms of four nested integrals, and the warp factor (</> 4 ) was expressed in 
terms of eight of them. 

Hence, the expressions presented in this paper correspond to the simplest solution to the 
space of SU{2) x SU{2) x Z 2 -invariant deformations of the Klebanov-Strassler background, 
and should include both the putative antibrane solution, as well as other perturbations of the 
Klebanov-Strassler field theory with various non-normalizable modes that have no R-charge (see 
for example [12, 13, 14]). 

In an upcoming paper we will use these results to perform the numerical integration of the 
putative anti-D3 brane solution of [4], allowing for the extra infrared singularity. This will relate 
the ultraviolet and the infrared integration constants and allow us to determine whether the ratio 
between the anti-D3 charge of the solution and the force on a probe D3 brane is the correct one. 

Note: While we were preparing this paper for submission we became aware of [15], in which 
some of the nested integrals obtained in [4] were explicitly performed. In the solution of [4] the 
most complicated field can be expressed in terms of 8 nested integrals, in [15] this is expressed in 
terms of 7 nested integrals, and in the present paper in terms of 3. Furthermore, as it is clear from 
the discussion above, we believe that it is premature to claim that a solution describes anti-D3 
branes without first computing the charge-to-force ratio and making sure it is the correct one. 
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2 Analytic Solution for Deformations of the Warped De- 
formed Conifold 

In this section we present our analytic formulae for the full set of 577(2) x SU (2) x Z2-invariant 
deformation space around the Klebanov-Strassler solution 4 . The method follows that of [4], but 
here we present numerous analytical improvements. 

We use the ansatz for the supergravity background fields proposed by Papadopoulos and 
Tseytlin (PT) [18], which is the most general ansatz consistent with the SU(2) x £77(2) x Z2- 
symmetry of the Klebanov-Strassler background: 

del = e 2 ^- dsl + e" 6 ^ dr 2 + {g\ + gl) + (gl + gl) + e" 6 ^ g\ , (1) 

where all the functions depend on the variable r. The fluxes and dilaton are 

H3 = 5 (k - f) g 5 A (gi Ag 3 + g 2 A g 4 ) + dr A (f g x A g 2 + k' g 3 A g A ) , 

F 3 = Fg 1 Ag 2 Ag 5 + (2P-F)g 3 Ag i Ag 5 + F'dTA(g 1 Ag 3 + g 2 Ag 4 ), (2) 

F 5 = J- 5 + *-F 5 , J r 6 = [kF + f(2P-F)]g 1 Ag 2 Ag 3 Ag4Ag 6 , 

$ = $(r), C = 0, (3) 

where P is a constant while /, k and F are functions of r and a prime denotes a derivative with 
respect to r. 

2.1 The Borokhov— Gubser method and the zeroth— order background 

The method introduced in [5] relies on the existence of a superpotential W whose square gives 
the potential, namely 

1 , dW dW 

v ^s G " b ww (4) 

The fields <j) a (a = 1, ...,n) are expanded around their supersymmetric background value in the 
form 

r = <j> a o + <K(X) + 0(X 2 ) , (5) 

where X represents the set of perturbation parameters, 0" is linear in them, and 0q are the 
functions in the Klebanov-Strassler solution, written explicitly below. The method amounts to 
splitting n second-order equations (n = 8 for us), into In first-order ones, out of which n of 
them (those for the conjugate momenta £ a ) form a closed set. The conjugate momenta £ a are 
defined as 

€. , GMo) (§ - *M *f) • m\ , \ » (g 6c ™) . (6) 



See [16, 17] for a rigorous derivation of a consistent truncation on T ' which includes several additional 
modes. 
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They are linear in the expansion parameters X and non-zero for non-supersymmetric solutions 
only The set (£ a , a ) satisfies the equations: 



^ + 6>M\(0 O ) 
d(Pl ~ M%(0 O ) <P\ 



dr 







i, ■ 



(7) 
(8) 



Note that equations (8) are just a rephrasing of the definition of the £ a in (6), while the equations 
in (7) imply the equations of motion [5]. The functions £ a should additionally satisfy the zero- 
energy condition 

4.^=0. (9) 
dr 

We will denote the set of functions <f) a , a = 1, 8 of the PT ansatz in the following order 



b a = (x,y,p,AJ,k,F,<S>). 



The field-space metric in (6) is 



G ab <p' a <p' b = e ip+iA 



2 4 



+-e- 



-*-2x ( e -2y f ,2 + e 2y k ,2 + 2e 2$ F ,2j 



and the superpotential is given by 
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W{<f)) = e AA - 2p - 2x + e 4A+4p coshy + - e * A+ * p ~ 2x (f(2P-F) + kF) 

2 

The background fields are given by the Klebanov-Strassler solution [1] 



(10) 



;i2) 



e x ° 

e 6po 

fo 

k 

Fo 
$ 



^/i(r) 1/2 (| sinh(2r) - r) 
tanh(r/2) , 

(| sinh(2r)-r) V3 



1/3 



24 



^i(r) sinh 2 r 



/i(r) (| sinh(2r) -r) 2/3 sinh 2 r 

(r cothr — 1) (coshr — 1) 
sinh r 

(r cothr — 1) (coshr + 1) 



3-2 9 
-P 



-P 



sinh r 



(sinh r — r) 
sinh r 



(13) 



0. 



The warp factor h cannot be found analytically. Its expression as an integral is given below in 
equation (23). 
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2.2 The first— order deformation: £ a equations 

As in [4] we shift to a slightly more convenient basis £ a , defined as 

g = (36 - 6 + U, 6, -36 + 26 - £ 4 , -3d + 6 - 2£ 4 , 6 + 6, 6 - 6, 6, 6) • (14) 

The equations in the order in which we solve them, are 5 

g = e~ 2 * [2Pf - F (Jo - fc )] g (15) 

g = -e" 2 ^ [2P/ - F (/o " fco)] 6 (16) 

Is = -\Pe- 2xo ii (17) 

§ = -h-\e- 2xo (P-F )£i (18) 

g = - sinh(2 y ) g - cosh(2 y ) g + 1 e" 2 ^ (/ - fc ) 6 (19) 

o 

g = (Pe 2 y° - sinh(2y ) F ) 6 + (P e 2yo - cosh(2 2/0 ) F ) 6 + £ (/o " to) 6 (20) 

g = 3 e - 2 *o-6po | 3 + [ 5 e -2.o-6 Po _ e -2, (2 p fo _ Fq (/q _ ko) } ] £ (21) 

g = g cosh sinh j/ (2 g + g + | 4 



+2 



( Pe 2i/o _ cosh (2y )F ) g+ (Pe 2 * - sinh(2 y ) F ) g . (22) 



i(r) = / — - — 3 , (24) 



The key development we present here is to solve for all the 6 in terms of two simple integrals, 
one of which is the KS warp factor: 

h(r) = h -32P 2 [ T U C ° th U 2 - 1 (cosh u sinh u-u) 1/3 du , (23) 
J sinh u 

du 

(cosh m sinh u — u 

with /io = 18.2373 P 2 a numerical constant. 

In solving the system of £ equations, we make the following key observations. In the equations 
for 6 and 6> we note that 

e -*«o[2P/ -F (/ -fc )] = ^. (25) 

This implies 

i 1 = X 1 h(r), Z4 = -X 1 h(r)+X 4 . (26) 
To obtain 6 we use the relations 

Pe 2 *' - sinh(2y ) Fo = ~\ (fo + h)' , (27) 
Pe 2 y°-cosh(2y )F = ~ (/„ - k )' , (28) 



3 We have accounted for two misprints in the published version of [4]. 



which yields 



g = —(fo + ko)' &-l(/ -* b ) / & + i(/o-*b) Ir- (29) 



Integrating by parts and using (18) and (26), we get 

§ = ~\ ( (/o + *b) - i ( (Jo - fcb) fe)' - \*x h' . (30) 
This easily integrates to 

Is = -5 (Jo + fco) Is - \ (fo - feo) £~e - ~ *i /i(r) + X 8 . (31) 



For £3 we observe that 



e -2* -6 P0= 4 sinh 2 r _ 1 e 3 , fa (32 x 
3 sinh 2r-2r3f 3h ' 



where £3^ is the solution to the homogeneous equation, namely 

£ 3th = sinh 2 r - 2r. (33) 

As a result, 

g = 3e- 2l »- 6 «6 + [5e- 2 ^- 6po - e" 2 * (2P/ - F (/ - fc„))] 6 

= 3e- 2x °- 6po i 3 + X 1 (-J^ _//) . (34) 

Henceforth, 

| 3 = -|x 1 /i(r) + |x 1 |3 J h / du^- + X 3 ^,h. (35) 

The above observations allow us to write the full solution in terms of two simple integrals. We 
collect here the full solution 

6 = X l h(r), (36) 

£3 = — X\ h(r) P 2 Xi csch 2 r (sinh r cosh r — r) 4//3 

3 3 

128 

P 2 Xi (sinh r cosh r — r) j(r) + 2X 3 (cosh r sinh r — r) , (37) 
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£4 = -X 1 / i (r)+X 4 , (38) 

£5 = -^-X 1 j(r)+X 5 , (39) 
~ 1 , . , cosh r sinh r — r , . . . . , 

«■ = -Sj77 a « W 2itoi7^ A 'M- < 40 > 

~ coshr , , . — 3 + cosh 2 r + 2 r cothr , , . , 

& = — r-T2-A 6 (r) + j- r-r A 7 r , 41 

sinh r 4 sinh r 

~ ~ t cotn T — 1 ~ 

£ 8 = P (r cothr -1) cothr£ 5 -P — £ 6 - -X, h(r) + X 8 , (42) 

smh r 6 

7 



where 



1 1 X 

A 6 (r) = X 6 + - (-r + cothr - r coth 2 r) | 5 (r) + - /i(r) , 

16 

A 7 (r) = X 7 — csch 2 r £5(7") H -P^i cschV (coshr sinhr — r) 1 ^ 3 

3 

64 

+-PX 1 j(r). 



(43) 



(44) 



Finally, £2 can be obtained through the zero-energy condition (9) or just by direct integration 
of (22). The latter will introduce another integration constant, X2, that one could then determine 
as some combination of the other ones via the zero-energy condition. We find 



2 1 

£ 2 = — 77 X 3 r cosh r + - X A cosh r + P X e cschr (coth r — r csch 2 r) 
o o 

+ PX 5 cschr (l — 2 r cothr + r 2 csch 2 r) + X 2 sinh r 
+ - P X 7 (-2 r coth 3 r + csch 2 r + r 2 cschV) sinh r 



108 



3 csch 3 r h(r) (6 r — 5 sinh 2 r + sinh 4 r) 



+ 2P 2 csch 5 r ( - 15 + 24 r 2 + 16 cosh 2 r - cosh 4 r - 32 r sinh 2 r + 4r sinh 4 r) 



x 



4 sinh 2 t j(r) — 6 (cosh r sinh r — r 



4/3 



(45) 



The zero-energy condition then amounts to 6 



Xo — — X 3 
3 



PX 5 --PX 7 = 0. 

5 2 c 



(46) 



2.3 The first— order deformation: (j) a equations 

The analytic expressions we obtain for the eight a modes are all double integrals, except for 4 
where we obtain a triple integral. This is a considerable improvement over all previous works. 
Depending on the reader's taste, the expressions may appear somewhat cumbersome but we find 
them crucial for explicit numerical computations which will appear in a companion paper [19]. 
To solve the system of <f) a , we also use a shifted basis [4] 

4> a = (x-2p- 5A, y,x + 3p,x-2p- 2A, f, k, F, $) . (47) 

6 The relation between the X a integration constants in this paper and the ones in [4] depends on the lower 
limit of the j(r) integral (24); it is not hard to check that this zero-energy condition is the same as that of [4]. 



The system of equations for the a modes is (in the order in which we actually solve them) 

& = -4e- 4A °- 4 «f 8 , (48) 

0' 2 = -coshy o 2 -2e- 4A °- 4 ^6, (49) 

0^ = _3e- 6 ^- 2 ^03-sinhy o 02-^e^ o ~ 4po (911 + 56 + 264) , (50) 

0; = 2e- 6 ^ 2a;o 03-sinhy o 2 + i e - 4Ao - 4 ^ (|i + 3| 4 ) , (51) 

e 2yo (F - 2 P) [2 2 + 8 ) + e 2w 7 - 2 e -4*>-4po+2*o+2 W ^ + |^ ? (52 ) 

>o (2 2 - 8 ) - 7 ] - 2 e - 4 ^-4 P o + 2xo- 2 ,o ^~ _ ; (53) 

" l(k-4> e + (k -f )fa)-2e- 4A °-^+ 2x ^ 7 , (54) 
^e- 2 ^[/ (2P-F ) + fc F ] (2^-203-504) +^e- 2 ^ (2P-F 



'5 



5' = e~ 2yo 



'7 



+ l -e^F ^+ l -e^ (A; o -/ o )0 7 -ie- 4A °- 4 ^6- (55) 

2.3.1 The 8 solution 

By directly integrating (48) and a little bit of massaging, we arrive at 
4> S =Y 8 -64X 8 j(t) + ^-!i(t) 

nA n ^ r ( U COthM - 1) , 

-64PX 6 / — ) —^du 

J sinh u (cosh u sinh u — u) 

2 ~ 16 

H /i(r) £5(7") H X\ csch 2 r (cosh r sinh r — r) 1 ^ 3 /i(t) 

P 3 

64 , N , N 64 r (sinh 2 u + l- u coth it) , s 

+ -X 1 h(r)j(r) + -X 1 / / ^7sM«)d«. (56) 

J J J sinh m (cosh u smh u — u) 

2.3.2 The 2 solution 

The solution to equation (49) is given by 

02 = cschr A 2 (r) , (57) 
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where 

/ T (—2 u coth 3 u + csch 2 u + u 2 csch 4 u) sinh 2 u 
~ 2/3^ du 
(cosh u sinh u — u) 1 

coth u — u csch 2 M , „„ „ f T 1 — 2 u coth u + u 2 csch 2 u 



-32PX 6 / ^du-32PX 5 / 273— du 

(cosh it sinh u — u) J (cosh it sinh u — u) 



32 /" T coshu sinhu 64 f T u coshu sinhu 

3 7 (cosh u sinh u — u) ' 3 J (cosh m sinh u — u) ' 

An 1/3 8 AA /" r 6u - 5 sinh2u + sinh4u , 

48 X 2 (cosh r sinh r - r) 1/3 + - X x \ - — ■ — - h(u) du 



9 J sinh 2 u (cosh u sinh u — m) 2//3 
32 2 _ /" r -15 + 24 u 2 + 16 cosh 2u- cosh 4 it - 32 u sinh 2m + 4m sinh 4 u 



-Fli / m du 

9 J sinh m (cosh u sinh u — u) 

+ 64 P 2 X /" r -15 + 24m 2 + 16 cosh2w - cosh4w - 32m sinh2w + 4m sinh4w 

27 j sinh 2 m (cosh m sinh m — m) 2 ^ 3 



(5J 



2.3.3 The 3 solution 

Equation (50) is solved by 



k(r) = — — ^L_-A 3 (r), (59) 
smh 2 r — 2 r 



where A 3 is specified as 

32 f T 112 /" T 
A 3 =y 3 — 3~"^ 4 / (cosh m sinh m — m) 1 ^ 3 c/m / ( cos ^- u sinhM — m) 1 ^ 3 h(u) du 

80 f T ~ f T 

— — (coshM sinhM — m) 1 ^ 3 £ 3 (m) du + 2 r cothr A 2 (r) — 2 / m cothM A' 2 (m) c/m . 



(60) 
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Expanding and simplifying this expression, it becomes 

/ T u cosh u (—2 u coth 3 u + csch 2 u + u 2 csch 4 u) sinh u 
; ^73 du 
(coshw sinhw — u) 

/ T u coth u ( coth u — u csch 2 w) 
273-^ du 
(cosh u sinh u — u) 

/ T u coth u ( 1 — 2 u coth u + u 2 csch 2 w) 
m du 
(cosh u sinh u — u) 

32 „ J _ [ T u cosh 2 u , f T .!/ 3 , 

+ — A 4 < 2 / — - du — I (cosh it sinh u — u) du 

3 I 7 (cosh it sinh u — u) ' J 



^- X 3 ^5 y (coshw sinhw — m) 4 ^ 3 + 4 J 



u 2 cosh 2 it 

du 



(cosh m sinh u — u 



2/3 



w coshw sinhw 
+ 2r cothr A 2 (r) + 64 X 2 / ^ du 



(cosh it sinhw — uf^ 

64 [ T 1/3 
+ — Xi / (cosh u sinh it — u) h{u)du 

H — — P 2 Xi y (cosh u sinh it — u) 4 ^ 3 j(u) du 

P 2 Xi I csch 2 u (coshu sinhu — uf^ du 



9 

,2„ 



16 r u coth u csch it (6 u — 5 sinh 2 u + sinh Au) 

— X 1 / 573 h(u) du 

9 J (cosh u sinh u — u) 

64 2 x „ /" T it cothu ( — 15 + 24 u 2 + 16 cosh 2 u — cosh 4 it — 32m sinh 2 u + 4u sinh 4m 



P Xl / i TTs dw 

9 J sinh m (cosh u sinh u — u) 

128 /" r 

-— P 2 Xi / (-15 + 24m 2 + 16 cosh2w - cosh4w - 32m sinh2w + Au sinh4w) 
u coth w csch 2 w 

X --T TT ( 61 ) 

(cosh u sinh u — u) 1 
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2.3.4 The 1 solution 

Next comes 0i which we express concisely in terms of A2 and 03 : 

40 2 160 f T 
0i =Yi + — X 4 j(t) - - 03 (r) — X 3 / (cosh u sinhu - u) 1/3 

+ y cothu A 2 (u) (in — ~ cothr A 2 (r) H — P 2 Xi J csch 2 u (coshu sinhu — u) 2 ^ 3 du 

10240 „ 2v f T , , . , a/3 v , 80 T /i(u) 

H — — P ai / (cosh w smh u — u) j{u)du——Xi / ^ du . 

81 7 27 7 (cosh m sinh m — w) 

(62) 

2.3.5 The (0 5 , 6 , 7 ) solutions 

The fields 05,6,7 are determined by a system of coupled ordinary differential equations. The 
homogeneous solutions are easily found and then we apply the Lagrange method of variation of 
parameters to find the following expressions 

1 1 

■ sech 2 (r/2) [r + 2 r cosh r - (2 + cosh r) sinh r] A 5 (r) + — A 6 (r) + A 7 (r) , 



(63) 

A 5 (r) + T -^A 6 (r)+A 7 (r), (64) 



2 1 + cosh t 

6 = r (2 3— — ) - coth(r/2) + sinhr 

\ 1 — cosh r J 

07 = (— coshr + r cschr) A 5 (r) — cschr A 6 (r) , (65) 
where 

1 „ / -, \ , 9 7 / n ^ f T ( u cothn — 1) csch 2 u ~ . , , 
A 5 =F 5 --P (r cothr - 1) csch 2 r0 8 (r) - 32 P / ^ —^( u )du 

7 (cosh u sinh u — u) 

1 f T [ T 2 + cosh 2 u 

+ — X7 / csch h [2 m (2 + cosh 2 u) — 3 sinh 2 u] h(u) du — Xq / 7 /i(w) 

4 7 7 sinh m 

/" T 2 r / 2 m 7 / \ z / \ In cosh r sinh r — r 

+ / csch it — 3 coth u + u ( 2 + 3 csch it) I h(u) ^(u) du P 7 A 2 (t) 

' 2 sinh 



1 Z" 7 " X\ f T 

— P / csch 4 u (cosh u sinh it — u) A 2 (u) — ^-7^ / (2 + cosh2 u) csch 4 u h 2 {u) du 

16 r 

— P Xi J csch 4 w [2 m (2 + cosh 2 m) — 3 sinh 2 m] j(u) h(u) du 
4 f T 

— PXi / csch 6 w (coshw sinhw — u) 1 ^ 3 [2u (2 + cosh 2 w) — 3 sinh 2 w] h{u) du , (66) 
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Y 6 - -P [-t + cothr + r (-2 + r cothr) cschV] 8 (r) 

'" r \—u + coth it + it (—2 + u coth it) csch 2 itl ~ , , 
-32P / ^ i ^r-L 



(cosh it sinh it — it 



2/3 



+ - X 7 J [cosh2it + csch 2 it (3 + 2u 2 — 6it cothit + 3 it 2 csch 2 it)] h(u) du 
+ Xq / csch 2 it [3 coth it — u (2 + 3 csch 2 it)] h(u) du 



+ J [l + (3 + 2 u 2 — 6 it coth it) csch 2 it + 3 it 2 csch 4 it] fa(it) dit 

- ~ P [2 coth 2 r (-1 + r cothr) + csch 2 r - r 2 csch 4 r] A 2 (r) 
1 /" r 

+ - P I [2 coth it (—1 + u cothit) + csch it — u 2 csch it] A 2 (it) rfu 



T f csch 4 u [—2 it (2 + cosh 2 it) + 3 sinh 2 it] 



12P 



1 « 

h(u) H Pcsctrit 

36 



x 



8j(u) sinh 2 it + 6 (cosh it sinh it — it) 1 / 3 — 28 + 32m 2 + (31 + 16m 2 ) cosh 2 1 



— 4 cosh 4 it + cosh 6 it — 48 it sinh 2 it 



h{u) du 



A 7 =Y 7 + P [-r + coth t + r (-2 + r coth r) cschV] </> 8 (r 
[—it + cothit + it (—2 + it cothit) csch 2 it] ~ 



64 P 



^s(m) du 



(cosh it sinh it — it) 2 ^ 3 
+ X 7 J [— 1 + (—3 — 2 u 2 + 6 it cothit) csch 2 it — 3 u 2 csch 4 it] h(u) du 

+ Xq / csch it [2 u (2 + cosh 2 it) — 3 sinh 2 it] h(u) du 



[— 2 — 2csch 2 it (3 + 2 it 2 — 6 it coth it + 3 u 2 csch 2 it)] h(u)^(u)du 
PcschV (l - 2r coth r + r 2 cschV) A 2 (r) 
P y csch 2 it (l — 2 it coth it + it 2 csch 2 n) h! 2 (u) du 
csch 4 it [2 it (2 + cosh 2 it) — 3 sinh 2 it] 



6P 



h(u) P csch it 

V ; 9 



8j(u) sinh 2 it + 6 (cosh it sinh it — u) 1 ^ 3 

[-9 + 16it 2 + 8 (l + it 2 ) cosh2it + cosh4it - 24 it sinh2it] J h(u)du. 
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2.3.6 The </> 4 solution 



While all the <f) a modes so far have been double integrals, we obtain for 04 a triple integral 
expression 



!>4 t 



+ 16 P 



3 



h{u) 



(cosh u sinh u — u 



,2/3 



du + 32 P 



r (w cothw — 1) csch uKr(u) 



(cosh u sinh u — u 



,2/3 



A 7 («) 



(cosh u sinh w — w 



2/3 



X 



5A 5 (u) + 2P -0 1 (n) + 3 («) 



32 
du 



P J (u cothu — 1) csch 2 u (coshu sinhu — u] 



1/3 



(69) 



It may be possible that this expression can also be reduced to double integrals, but we could not 
find any obvious way to do it. This completes the solution to the system. 



3 Boundary Conditions and Anti— D3 Branes 

The deformation space we have solved for is a fourteen-dimensional linear space 7 and has nu- 
merous solutions, out of which one has to fish out the possible solution for backreacted anti-D3 
branes by imposing appropriate boundary conditions. 

The strategy for doing this was explained in detail in [4]: one should eliminate all integration 
constants that give divergent fields in the IR or non-normalizable modes in the UV. Furthermore, 
to argue that the solution corresponds to D3 branes one should set the divergence in the warp 
factor perturbation (given by ^4) to be commensurate with the divergence coming from the five- 
form. This should fix all the integration constants in terms of the would-be anti-D3 charge. 
Nevertheless, to obtain the precise values of these constants, one needs to relate the infrared and 
ultraviolet expansion parameters of the modes presented in this paper, and this can only be done 
by numerical integration. 

In particular, this calculation will fix the constant Xi in terms of the additional anti-D3 
brane charge iVgg which is given by the UV value of the mode 5 + <p & . Having done that, one 
can test whether the candidate solution describes indeed anti-D3 branes by using the universal 
form of the force exerted on a probe D3 brane by any solution to the equations [4, 11]: 

P D3 = - Xl[Nm] (70) 

3 (sinh rcosh r — r) 2//3 

and comparing this with the force anti-D3 branes should exert [3, 11]. The expressions for 
the 4> a ) obtained here in terms of single and double integrals turn out to be crucial in this 
endeavor, on which we intend to report soon. 
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7 One of the X's can be eliminated through the zero-energy condition (46), while Y\ is a gauge degree of 
freedom corresponding to a constant rescaling of the 4D metric. 
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